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Explosive synchronization can be observed in scale- free networks when Kuramoto oscillators have 
natural frequencies equal to their number of connections. In the current work, we took into account 
mean-field approximations to determine the critical coupling of such explosive synchronization. The 
obtained equation for the critical coupling has an inverse dependence with the network average 
degree. This expression differs from that calculated when the frequency distributions are unimodal 
and even. In this case, the critical coupling depends on the ratio between the first and second 
statistical moments of the degree distribution. We also conducted numerical simulations to verify 
our analytical results. 
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INTRODUCTION 

Synchronization of coupled oscillators has been in- 
tensively studied because of its ubiquity in the real 
world [TJ [2J. When a collection of oscillators are cou- 
pled as a network, it can be observed the emergence of a 
synchronous state (21 E]. Such onset of coherent collective 
behavior has been verified between neurons in the central 
nervous system, communication networks, power grids, 
social interactions, animal behavior, ecosystems and cir- 
cadian rhythm [2J. 

The level of synchronization of a system is the conse- 
quence of a combination of the type of oscillators, the 
connectivity organization, the time-delay and the inter- 
action function |TJ [2J . Particularly, the network topology 
has a strong influence on the value of the critical cou- 
pling [4-8 and on the stability of the fully synchronized 
state [21I5HTT]. For instance, Watts and Strogatz [T2] ver- 
ified that the decrease in the average shortest path length 
in small-world networks facilitates a more efficient cou- 
pling and therefore enhances the level of synchronization. 
In addition, Nishikawa et al. [TT] suggested that networks 
with an homogeneous degree distribution are more syn- 
chronizable than heterogenous ones. 

The network structure is not only important to en- 
hance the level of synchronization, but also to permit 
the occurrence of phase transitions. Indeed, many works 
have verified second-order phase transitions in networks 
of Kuramoto oscillators [2J . Recently, Gardehes et al. |13j 
showed that a first-order nonequilibrium synchroniza- 
tion transition can occur in scale- free networks. They 
suggested that this event is a consequence of a positive 
correlation between the heterogeneity of the connections 
and the natural frequencies of the oscillators [T3] . First- 



order phase transitions were also obtained experimen- 
tally and numerically by considering a network of Ros- 
ller units [TJ. Indeed, such phenomena is attracting the 
interest of many complex networks researchers {e.g. |14f - 

W)- 

Although the explosive synchronization has been ob- 
served in scale- free networks, the analytical expression 
that describes the critical coupling has not been deter- 
mined yet. In the current work, we obtained such ex- 
pression by considering mean-field approximations. We 
verified that the critical coupling has a inverse depen- 
dence with the network average degree. Our analytical 
results are compared with numerical simulations. 

The Kuramoto model considers a set of N oscillators 
coupled by the sine of their phase differences and phase 
oscillators at arbitrary frequencies [T7] . Each oscillator is 
characterized by its phase 9i(t), i = 1, N. In complex 
networks, each oscillator i obeys an equation of motion 
defined as 

dO N 

^=LJ i + Xj2 A ijM^-^i), i=h---,N, (1) 

i=l 

where A is the coupling strength, u>^ is the natural fre- 
quency of oscillator i, and Ay are the elements of the 
adjacency matrix A. so that Aij = 1 when nodes i and 
j are connected while — otherwise. The general 
Kuramoto model considers a random distribution of the 
natural frequencies and phases according to a specific dis- 
tribution g(uj) [TJ[2]. In most of the cases, the frequency 
distributions are unimodal and symmetric around a mean 
value wo 0. 

In the current work, we considered a modified ver- 
sion of the Kuramoto model as proposed by Gardehes 
et. al. |13j . More specifically, the natural frequency u>i 
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of the node i was assigned to be equal to the node's de- 
gree ki, i.e., uii = ki. Therefore, g(uj) = P{k), in which 
P(k) is the degree distribution of the number of connec- 
tions. This choice for the frequency distribution leads to 
the explosive synchronization in scale-free networks |13| . 
Gardehes et. al. verified that this effect is due exclusively 
to the positive correlation between the network struc- 
ture and the dynamics. When this correlation is broken, 
a first-order transition is no longer observed, whereas a 
second-order transition occurs |13j . 

In order to analyze the interplay between structure and 
dynamics in the Gardehes et. al. model, we considered 
the mean-field approach proposed by Ichinomiya |18j . 
First, we characterized the network by its degree dis- 
tribution P(k) and introduced the density of the nodes 
with phase 9 at time t for a given degree k, denoted by 
p{k;9,t), which is normalized according to 



p{k;9,t)d9 = 1. 



(2) 



The continuum limit of Eq. [T]is taken by considering the 
absence of degree correlation between the nodes in the 
network. Observe that this is a typical assumption in 
mean-field approximation pQ. In this regime, the proba- 
bility that a random edge is attached with a node with 
degree k and phase 9 at time t is given as 



kP(k)p(k;9,t) 
(k) 



(3) 



where (k) is the network average degree. Replacing 
U!i = ki in Eq. [T] and taking the continuum limit in 
the mean-field approach using Eq. [3] we obtained 

d0(<) -k + Xk f dk' [ d6' k ' P S^ p(k; 6', t) sm{6 - 6'). 



dt 



(k) 



(4) 

The order parameter, which quantifies the level synchro- 
nization of the network, is defined as [T5I [T§] 

re^ (t) = J dk J d9kP{k)p(k;9;t)e ie , (5) 

where < r < 1 and tp(t) stands the average frequency 
of the oscillators. 

Multiplying EqjH] by e zS , taking the imaginary part 
and including in Eq. [4j we obtained 



d9 

— = k + XkrsinOib — 9), 
dt 



(6) 



which is the Eq. [4] written in terms of the order parame- 
ter. 

In order to let the equations of motion in function of 
known parameters of the network, we set a reference ro- 
tating frame ip(t) — fit, where Q is the average frequency 
of the network. In the case of the Gardehes et al. model. 



i.e. g(ui) — P(k), the average frequency is equal to the 
network average degree (Q = (k)) |13| . Defining a new 
variable as <f>(t) = 9{t) — ip(t) and replacing in Eq. [6j we 
obtained 



dt 



= (1 — Xr sin <f>)k — (k) . 



(7) 



We redefined the density of oscillators p in terms of 
the new variable 4>, i.e. p = p{k;<j),t). This density of 
oscillators must satisfy the continuity equation |18j 

dp(k;<j,,t) d 

— dt — gilw;^*)} = °> ( 8 ) 

where v$ = Since we were interested in the analysis 
of the steady state of the system, we obtained the time- 
independent solutions of Eq. [8] i.e. 



p{k\(j>) 



8 [4> — arcs 

A(k) 

. \(k-(k))- Xkr sin <j>\ 



if 



\k-{k)\ 



< Xr, 



otherwise, 



(9) 

where §(■) is the Dirac delta function and A(k) is the 
normalization factor. The first solution is respective to 
the synchronous state, i.e. t| = 0, corresponding to the 
oscillators which are entrained by the mean field. On 
the other hand, the second one is the density of the non- 
entrained oscillators, i.e. p{k;4>) ~ j^-j j3J[TB]. Thus, to 
compute the integrals in Eq. [5] we redefined it in terms 
of the variable 4> and separated the contribution of en- 
trained and non-entrained oscillators 



dk 



|fc-WI 



dk 



fc-(fc) 



>\r 



P(k)kp{k; ^)e**#. (10) 



Rewriting the second integral in Eq. 10 and noting that 
p(k; <ft) is 7r-periodic in (f>, we obtained 



o J(k)/(i-\r) 2ir(k- (fe) - kXramfi) 

„ A • 27r((fc> -fc + fcArsin0) Y 



where k mm is the minimum degree in the network. Thus, 
only the contribution of the oscillators entrained in the 
mean-field is accounted in the summation of Eq. |10| 



(k)r 



(k)/(l-\r) 



exp 



'(fe>/(l+Ar; 
x kP(k)dk 



1 

Ar 



(k) 



From the imaginary parte of Eq. we obtained 

'k- (k) 



(fc>/(l-Ar) i 
kP(k) 

(k)/(l + \r) 



Xr 



dk = 0, 



(11) 



(12) 
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and from the real part, 



,(fc)/(i-Ar) / 1 /h-lkW 2 

{k)r= fc p(A.)Ji_ J_(^W\ dk _ 

J(k)/(1+Xr) V A r \ k J 

(13) 



Considering a; = (fc — (fc))/Ar, we obtained 

,(fe>/(l-Ar) 

(fc)r = Ar / P(Ara;+ (fc))(Arx+ (fc)) 

W/0:Ar) 



X 4/1- 



^Arx + (k) 
For r ^ and letting r-)0 + ; 

(k) = xj { %((k))(k) ] ji~- 

we reached to the critical coupling 



dx. 



(14) 



n(k)P((k)Y 



(15) 



(16) 



Therefore, the critical coupling presents an inverse de- 
pendence with the average network degree and P((fc)). 
This dependence is very different from that observed 
when it is taken into account other types of frequency dis- 
tribution g(ui). For instance, if g(ui) is symmetric about 
a single local maximum ujq (e.g. ujq = 0) the critical 
coupling is given as [151 [TU] 



Al°) = 



(k) 



7T 5 (0) (fc 2 ) 



(17) 



Thus, for scale- free networks, as N — > oo the critical 
coupling Ac become smaller, since the ratio (k) / (k 2 ) 
diverges. On the other hand, for g(uj) = P(k), this effect 
should not be observed when N — > oo, because the crit- 
ical coupling depends only on the average degree (k) of 
the network. 

In order to check the validity of Eq. [9] and [TBJ we 
considered numerical simulation. We increased the cou- 
pling strength A adiabatically and computed the sta- 
tionary value of the global coherence r for each value 
Ao, A + SX, Ao + nSX, with increments SX = 0.02, as 
done in [13] ■ Fig. [T] shows the dispersion of the phases <f> 
as function of the node's degree k for a BA network with 
N = 10 3 nodes with (k) = 6. As we can see in this figure, 
for A = 2.0 the system starts to present partial synchro- 
nization, suggesting that the critical is between A = 1.0 
and A = 2.0. Note that for A = 4.0, the numerical results 
of the phases <j> are in good agreement with the theoretical 
solution, specially for the highly connected nodes. Fig. [2] 
also presents the dependence of the phases <f> on the de- 
gree k for a Erdos-Renyi (ER) network with N = 10 3 and 
(k) = 6. As in Fig.[l] we observed the same behavior for 
the ER network, as the coupling A becomes higher, the 
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FIG. 1: Distribution of the phases cj> as a function of the 
degree k for a BA network with iV = 3 • 10 3 nodes and average 
degree (k) = 6. 



phases approaches the theoretical solution. Therefore, 
our results suggest that the solution of p(k;4>), given by 
Eq. [9j is valid. 

Once we have verified the validity of Eq. [9] , we esti- 
mated the critical coupling considering numerical data. 
We took into account an ensemble of N net networks, 
{Ni, N2, N net }, with the same number of nodes iV" 
and same average degree (k) = 6. We estimated the 
critical c oupl ing A c as an average over this ensemble by 



using Eq. 

. (JV=2000) 
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In this way, we obtained A 

(Af=3000) 



(jV=1000) 



0.65, 



Ar ; = 0.65 and A^ 1 " ' = 0.66 for the ER net- 
works. Fig. [3] shows the coherence diagram of r as func- 
tion of A for ER networks with N = 10 3 , 2 • 10 3 and 3 • 10 3 
nodes. As we can see in this figure, the critical coupling 
A c does not depend extensively on the total number of 
nodes N in the network, since the theoretical estimation 
for the critical coupling in Eq. [16] depends only on the 
network average degree. 

In order to conducted numerical simulations to ver- 
ify the validity of Eq. |16| we considered the model by 
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FIG. 2: Distribution of the phases <j> as a function of the 
degree k for a ER network with N = 3- 10 3 nodes and average 
degree (k) = 6 . 
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FIG. 3: Synchronization diagram for ER networks with for- 
ward continuation of the coupling strength A with steps of 
S\ = 0.02. The networks have N = 1 • 10 3 , 2 • 10 3 and 3 • 10 3 
and the same average degree (k) —6. The black bar depicts 
the theoretical value for the critical coupling A c . 
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FIG. 4: Synchronization diagram for BA networks with for- 
ward continuation of the coupling strength A with steps of 
S\ = 0.02. The networks have N = 1 • 10 3 , 2 • 10 3 and 3 ■ 10 3 
and the same average degree (k) —6. The black bar depicts 
the theoretical value for the critical coupling A c . 
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FIG. 5: Coherence diagrams for configurational models with 
degree distribution P(k) ~ fc -7 : (a) 7 = 2.4 (b) 7 = 2.5 (c) 
7 = 2.6 and (d) 7 = 2.7, with forward continuation of the 
coupling strength A with steps of SX = 0.02. 



Barabasi and Albert (BA) and a configuration model. 
Networks generated by the BA model are characterized 
by a distribution of connections following a power law, 
i.e. P(k) ~ fc -3 [20]. The configuration model allows 
to generate networks with a given degree sequence |21| . 
Fig. [4] shows the coherence diagram for a BA network. 
Performing the same procedure described above to es- 
timated the critical coupling, we obtained the values 

A ( M 0) ^ Lgl) A (JV=2000) ^ 1M ^ A ( W =3000) ^ 

1.52, which are in good agreement with the results from 
the numerical simulation. Also, Fig. [5] shows the explo- 
sive synchronization in scale-free networks with degree 
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distribution P(k) ~ fc~ 7 constructed using the configu- 
rational model [2J with 7 = 2.4, 2.5, 2.6, and 2.7, con- 
sidering degrees in the range 1 < k < 100. 

When the frequency distribution, g(uj), is unimodal 
and even, the critical coupling tends to vanish as TV — > 00. 
On the other hand, the assumption that the natural fre- 
quencies are equal to the degree implies that the critical 
coupling does not suffer significant variations. This fact 
can be observed in Fig. |4j In addition, we did not ver- 
ify that for the forward continuations of A, the critical 
coupling increases extensively with the number of nodes. 
This result was obtained in [T3], where authors consid- 
ered a star network as an approximation of scale-free net- 
works. Therefore, although star networks exhibit the first 
order phase transition, the critical coupling does not have 
the same behavior as verified in scale-free networks. 

The analysis performed in the current work helps to 
understand the relationship between the structure and 
the explosive synchronization in scale-free networks. The 
obtained expression for the critical coupling does not de- 
pends on the ratio (fc) / (fc 2 ), as observed in the case 
when g(uj) is symmetric. Indeed, the obtained critical 
coupling has a inverse dependence with the network av- 
erage degree, (fc), and P((fc)). 
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